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Braneworld scenarios with compact extra-dimensions need the volume of the extra space to be 
stabilized. Goldberger and Wise have introduced a simple mechanism, based on the presence of a 
bulk scalar field, able to stabilize the radius of the Randall-Sundrum model. Here, we transpose 
the same mechanism to generic single-brane and two-brane models, with one extra dimension and 
arbitrary scalar potentials in the bulk and on the branes. The single-brane construction turns out 
to be always unstable, independently of the bulk and brane potentials. In the case of two branes, 
we derive some generic criteria ensuring the stabilization or destabilization of the system. 
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Models with compact extra-dimensions need the vol- 
ume of the extra space to be stabilized. This is a typical 
problem of string-theoretical constructions, in which the 
existence of a mass term for moduli is vital for the via- 
bility of the models. The inevitability of 9 or 10 spatial 
dimensions in string theory provides one of the main mo- 
tivations for the recent interest about extra-dimensions, 
although many scenarios that are not necessarily linked 
to string theory have also been envisaged. Among these, 
the one proposed by Randall and Sundrum Q has been 
object of intense study. One of the (many) features that 
make the scenario so interesting is the presence of a single 
extra-dimension, that allows for a detailed - and not too 
complicated - analysis of the model and of its possible 
modifications. 

Soon after the proposal of 0] (that was not address- 
ing the problem of the stabilization of the extra-space), 
Goldberger and Wise Q have considered a simple mech- 
anism, based on the presence of a bulk scalar field, able 
to stabilize the radius of the Randall-Sundrum model, 
as was proved in detail in and 0. Since then, the 
literature has been enriched by many works about brane 
models in presence of bulk scalars. The importance of 
radion stabilization is related to the phenomenological 
consistency of the models @ as well as to possible obser- 
vational consequences 0. For some specific form of the 
scalar field potential the behavior of the system could be 
reconstructed Q, by solving the bulk Einstein equations. 
In most cases, however, such solution can be obtained 
only cither in the limit in which the backreaction of the 
scalar field on the background geometry is negligible, or 
numerically. For instance, the dynamics which leads var- 
ious examples of two-brane models to stable equilibrium 
was recently studied numerically in Q . 

In the present paper we will analyze the spectrum of 
scalar perturbations for (4 + l)-dimcnsional brane mod- 
els in the presence of a bulk scalar field. Our main goal is 
to provide criteria that can tell us in which case tachyon 
modes emerge: a system with tachyon modes is unsta- 
ble and would quickly decay into some different ground 
state. We analyze the equations describing scalar pertur- 
bations around a background with (3 + 1) -dimensional 



flat branes, without referring to any specific potential for 
the scalar field, and without making any approximation. 
We classify models only on the basis of the behavior of 
the background metric and scalar field. We find that, 
under general assumptions, models with only one brane 
always feature tachyon modes, consistently with the con- 
jecture in 0| . In the case of models that include two 
branes, we provide general criteria about the presence 
or the absence of tachyon modes in the spectrum of the 
theory. 



I. SETUP AND MAIN EQUATIONS 

We consider a five-dimensional spacetime, with four- 
dimensional Minkowski sections (spanned by the coor- 
dinates x^) and one compact extra-dimension y. Two 
three-branes are located at the fixed points y = j/_ 
and y = y+ > y~ of a Z2 symmetry y «-» 2 y — y + , 
y «-> 2 y_ — y. In section II we will specialize to the case 
in which only one brane is present. The matter content 
of the model is given by a bulk scalar field tp, with a bulk 
potential V (p) and additional potential terms U± (if) lo- 
calized on the branes. A possible cosmological constant 
in the bulk and eventual brane tensions are included in 
the definition of V and U± . The (background) metric for 
such system can be written as 

ds 2 = 9mn dx M dx N = a (y) 2 (n a p dx a dx^ + dy 2 ) , 

(1) 

with r] a/3 = diag(-l, 1, 1, 1). 

Denoting by a prime the derivative with respect to the 
extra coordinate y, and defining H (y) = a' (y) /a (y), we 
write the background equations of motion as 



dV 



dU± 



<p'o+3Htp' -a — (ipo) = S(y±) a— — (ip ) , 



dip 



d(p 



H' = H 2 -^p' 2 ~S(y ± )^ K 2 U ± (p ), (2) 



with the constraint 

K 2 

6 



-^-a 2 V(p ) 



(3) 



2 



k 2 represents the five-dimensional Newton constant. 

The Israel junction conditions, together with the inte- 
gration of the Klein-Gordon equation across the branes, 
give the boundary conditions at the brane locations 



H(y±) = ±—a(y±) U±(<po{y±)) 



(<Po (y±)) 



(4) 



In order to study the stability of this system, we per- 
turb the metric and the scalar field as gMN — 9mn + 
h M N, V = <Po + Sip, with 



h a p = 2a(y) (tp r) a(i - dadpE) 
h a5 = h 5a = -a (y) 2 d a B 
h 55 = 2a(yf (/) . 



(5) 



Since we expect possible instabilities to come from (3 + 
l)-dimensional scalar modes (the radion), we limit our 
analysis to scalar perturbations. 

The system turns out to be fully described by the fol- 
lowing three gauge invariant quantities 

* = V + H (B — E') , 

$ = t/> + H (B - E') + (B - E')' , 

A ee Sip + tp' (B - E') , (6) 

and the linearized bulk Einstein equations reduce, in 
terms of these gauge invariants, to the constraints 



$ = -2*, 
-<p' A = 3 (#' + 2ff*) , 



(7) 
(8) 



and to the dynamical equation (written here only in the 
bulk, we will treat the boundaries separately) 

□ * + 4ff*' + (6H 2 + 2H')V = l[<p' A'- a 2 ^(<p ) A], 
v ' 3 dip 

(9) 

from which we derive a wave equation for ^ in the bulk 



(10) 



where the box denotes the four-dimensional D'Alcmbcrt 
operator. 

The last equation can be rewritten as 



Uu ■ 



u = 0, 



H 



7 3 / 2 in' 



(11) 



where we defined, in analogy with the usual four- 
dimensional cosmological case u = (a 3 / 2 /^) ^. 

In general, there can be at least one point y^ in which 
f'o (yo) = an( i the equation of motion for ^ and u are 
singular. However, it is possible to go across the singu- 
larity using a third quantity whose equations of motion 



remain regular |12|. This is the equivalent (for our 5- 
dimensional setup) of the Mukhanov variable 0], i- e - 



,3/2 



A 



H 



(12) 



that obeys 



v" + nv V = , Z EE 1/0 , (13) 

z 



where the ratio 



z 
z 
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19/^ 



(i 



/ 2 



d 2 V 
dip 2 



2k a p' 4 
9 H 2 



dV 



4« 2 <p'q 2 _ 
3 H a dip 



(14) 



is always regular as long as H ^ 0. Going back to the 
metric perturbation is straightforward using the relation 



□ \I> = 



a 3 / 2 3 



(15) 



The boundary conditions for ^ can be obtained by 
perturbing the Israel junction conditions and the Klein- 
Gordon equation. Besides the perturbed quantities de- 
fined above, one has in this case also to take into account 
the perturbation £± of the position of the branes. Defin- 
ing the gauge invariant Z± = (± - [B (y±) - E' (y±)], 
the junction conditions yield Z± = 0, whereas the Klein- 
Gordon equation gives the following boundary conditions 

9± (V±) + 2 H (y±) * (y±)] + (y±) = , (16) 

where we have defined 

g± = {v±) rt^yT) 2 a {y±) {vo {y±)) ■ 

(17) 

Fourier-transformation of the above equations along 
the four ordinary dimensions amounts to the substitution 
of the □ operator with the mass eigenvalue m 2 . Solv- 
ing eq. (|1(J[) with the boundary conditions i|16|) leads to 
a discrete mass spectrum, corresponding to the various 
Kaluza-Klcin cigenmodes. In the next sections we will 
study the conditions leading to the generation of tachyon 
modes (to 2 < 0) that signal an instability of the model. 



II. ONE BRANE 

It possible to construct models that include only one 
brane. In this case, the Z 2 symmetry emerges naturally, 
with all quantities remaining regular in the second fixed 
point. In this section, we will locate the brane at y = y+, 
and we will assume the other fixed point y_ to be at 
j/_ = 0. In this case, the Zi symmetry at y = imposes 



<p' (0) = if(0) = 0. 



(18) 
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The constraint equation then implies that the poten- 
tial V (ip) has to vanish at y = 0. Close to y = 0, the 
solution of the Einstein equations (0 reads 



vo (y) 
H{y) 



^(0) + ia(0) 2 |^ (0))y 2 + O(y 4 ) 



dip 



(Vo(0)) 



y 3 + O (y 5 ) . (19) 



We see from the last equation that H (y) has to be nega- 
tive close to y = 0. In fact, this turns out to be the case 
everywhere: 



V ye] ;2/+ ] 



H(y)<0. 



(20) 



Indeed, if H reaches zero at some y = yn' 

(i) either ip' (jjh) does not vanish, and then, by eq. Q, 

H' (y H ) < ; 

(ii) or ip' (jjh) vanishes (but then, tp' ' is non-zero, other- 
wise the bulk scalar field would reduce to a cosmological 
constant). In this case, the derivatives of eq. show 
that W (y H ) = H" (y H ) = 0, while 



H'" (y B ) 



|« 2 ^ ,2 (te)<o. 



So, in both cases (i) and (ii), H should be turning from 
positive to negative, which is in contradiction with the 
fact that H (y < y H ) < 0. 

We will prove the existence of at least one tachyon 
mode in the perturbations of this setup. To begin with, 
let us define a function F whose zeros will give the 
Kaluza-Klcin mass spectrum: 



where 



F (to 2 ) ee g+ U m 2 (y + ) + m 2 * m 2 (y+) , (21) 



II m2 (y) ee *' m2 (y) + 2H (y) 4< m 2 (y) . (22) 



In eqs. ^ m 2 (y) is a solution of eq. (fTUf) that 

satisfies the boundary condition (|16fl at y = 0. Such 
boundary condition, as well as the regularity of eq. i|10|) 
in y = 0, imposes that ^' m 2 (0) = for every value of 
m 2 (this can also be seen from eqs. 1181 1 81) )■ Without loss 
of generality we set \E' m 2 (0) = 1 [lf|. Also note that 
eq. HI (J| l can be conveniently re-expressed as a Hamilton- 
like system 



-H 



n„,2 -2ff* 



n m 2 + -m 2 + 



2k' 



/ 2 
^0 



(23) 



with the following initial conditions in y = 0: 

(* ro2 ,n m2 ) = (l,0), (^ n2 ,n^ 2 ) = (0,m 2 ). (24) 

It turns out that F(m 2 ) always goes through zero for 
some negative value of to 2 . To reach this result, one 



should first note that the equation of motion for \E' m 2 
can be explicitly solved for to 2 — ► — oo. In this limit, 
eq. I|13(l reduces to 



ii i 6 

V m 2 + TO V m 2 V m 2 







(25) 



We can solve this equation imposing the boundary con- 
dition ^ (0) = 1. The solution reads 



v m 2 =-3 



o(0) 



-l/2 m 2 



sinh (|to| y) 



1 - 



_ 2 2 

Tfi y 



cosh (\m\y) 



\m\y 



(l + 0(m- 2 )) . (26) 



From this expression we can obtain the behavior of ^ m 2 
for |to| y 3> 1, 



*m"(l/) 



1 o(0) 



3/2 



(y) 



n m2 (y) = -V 



2 a(y) 3 / 2 a(0) 2 dF/^(^o(0)) 

l-l)^(M-) 

m| a(0) 3 / 2 <p> (y) 



J m \y . 



2 a(y) 3 / 2 a(0) 2 dF/^((^o(0)) 

y 



3H 



0(| 



(27) 



In particular, for y = we have in the limit |m| y+ ^> 1 

H a(0) 3/2 <^o(y+) j m u + 



n„ 



2 a(y + f 2 a(0)W/d^o(0)) 
\m\ *,„2 . (28) 



In the above equations, dV/d(p (y>o(0)) has the same sign 
as (p' (y) for y close to y = 0. If we assume that ip' van- 
ishes in an odd number of points in the range ] 0, y+[, this 
sign will be opposite to the one of ip' (y+): so, ^ m 2 (y + ) 
and IT m 2 (y+) will be both positive. This implies that 



F(—oo) = lim m 2 \E' m 2 (y+) = — oo . 



(29) 



In the opposite case, i.e. when Lp' Q vanishes in an even 
number of points in the range ] 0, y + [, we get that 
F(— oo) = +oo (of course, in the case "even number of 
zeros" , we include the configurations with no zeros at 
all). 

The equation of motion i|10|) or the system (|23|l can be 
solved explicitly in a second case, when m 2 = 0: 



*o(y) = 1-2 



H(y) 
a(yf 



a (yf dy , 



no(y) = f V a(yfdy. 

a (y) Jo 



(30) 



This solution is such that ^ (y) > 1, (y) > and 
IIq (y) > 0. So, F(0) — g + U (y + ) is of the same sign as 



4 



n 



4D- 

\ 111! 


F(m 2 ) 








FIG. 

™2 



0,mf 



The phase-space trajectory for the three modes 
oo, for a particular single-brane background 



with V(<p) = -Vo + (m 2 /2)ip% U(<p) =U - (fi /2)tp, and 
¥>o(0) = (2Vb/m 2 ) 1 ' /2 . In this example, there are no points 
where tp' = 0. The zero-mode trajectory is the very short 
red line starting from (1,0) and moving upward and right. 
The mode m 2 is defined in the text as the first one hitting 
the axis = 0, when m 2 decreases from to — oo. On the 
upper right, we also show the function F(m 2 ) for the same 
model, in the case g+ > 0. The three black dots correspond to 
F(m 2 ) = 0, i.e., to the first three Kaluza-Klein modes. Two 
of them are tachyon modes. 



g + . On the phase-space diagram of Fig. ^ we plot the 
two solutions m 2 — > — oo and m 2 = obtained numeri- 
cally for a particular choice of background. 

For the last step of the proof, we should distinguish 
between three cases: 

(a) if tp' has an even number of zeros in the range ] 0, y+[ 
and c/+ is negative, or ip' has an odd number of zeros and 
g+ is positive, we immediately see that -F^O) and F(— oo) 
have opposite signs, so F(m 2 ) goes through zero for some 
negative value of m 2 0| , and the system has at least one 
tachyon mode. 

(b ) if ip' has an even number of zeros in the range ] 0, y+ [ 
and g+ is positive, or ip' has an odd number of zeros and 
g + is negative, F(0) and F(— oo) have the same sign, so 
we cannot draw an immediate conclusion. However, it is 
not difficult to show that there exists at least one value 
m\ < such that F(m\) and F(0) are of opposite sign. 

The crucial point is that in the limit m 2 — > — oo, the 
trajectory (^! m 2(y), Tl m 2(y)) always includes some points 
where >F m 2(y) is negative. When ip' has an even number 
of zeros, this is obvious, since we know that 
must be negative in the large — m 2 limit. When ip' has 



an odd number of zeros, we will focus on the first of them, 
i.e., the smallest coordinate value yo for which tp' (yo) = 
0. The sign of ^ m 2(y ) can be extracted in the limit 
of large — m 2 using eq. (|27fl . The leading nonvanishing 
term comes from O'/O, that diverges as y — > yo, but is 
multiplied by the overall factor ip' (y). We are thus left 



with 



(yo) 



■(0) 



3/2 



<p'o (yo) 



yo 



a(y f 2 a(0) 2 dV/d<p(M0)) 



(31) 



In this expression, dV/d<p(ipo(0)) has the same sign as 
ip (y) for y close to zero, whereas ip' ' (yo) has the opposite 
sign. As a consequence, >F m 2 (yo) turns out to be negative 
for large enough — m 2 . 

This fact has important consequences. When m 2 
varies from zero to — oo, the phase- space trajectory 
(^E' m 2 (y), II m 2 (y)) moves in a continuous way from solu- 
tion J3D|| to solution (O. The former is such that * > 1, 
while the latter includes points with 'F < 0. So, unavoid- 
ably, there exists a value m 2 for which the trajectory just 
touches the axis *F = for some value y: 



Vy, * m? (y)>0 and * TO? (y)=0. 



(32) 



In y = y, the derivative ^' m2 cannot be positive, since the 
trajectory comes from the region 'F > 0. So, H m 2(y) = 
^' m 2 (y) can be only negative or zero. Let us discuss these 
two cases: 

• if II m 2(y) < 0, we find that y must be equal to y+, 
because otherwise ' m 2 (y) would reach some negative val- 
ues for y > y, in contradiction with the definition of m 2 . 

• if H m 2(y) = 0, the trajectory reaches the origin of 
phase-space. Let us suppose that y is a special point 
where tp' vanishes. Then, the behavior around the ori- 
gin can be deduced directly from the perturbed Einstein 
equations JHJ, ©: it is easy to show that ^'^(y) van- 
ishes, but not ^'^(y). So, 

^! m 2 i (y)^a(y-y)\ IL^ (y) ~ 3 a (y - y) 2 , (33) 

where a is an unknown parameter. This behavior cannot 
occur for the mode m 2 , because , F m 2 (y) would change of 
sign in y = y, in contradiction with the definition of 
mf. So, y can only be in the range ]y n ,y+] where y n 
is the last point in which tp' vanishes. In this range, 
the trajectory is governed by the system which tells 
that the trajectory will remain in (^ m 2 , IT m 2 (y)) = (0, 0) 
until y = y+. 

So, we have proved that in any case, the mode is 
such that 



*m?(2/H 



0. 



n m 2( y+ ) <0. 



(34) 



The trajectory corresponding to m\ is shown in Fig. ^ 
for a particular choice of background. We finally obtain 
that F(m 2 ) = <7+IT m 2 (y+) is of sign opposite to <?+. Since 
we know that F(0) is of the sign of g + , there is always 
at least one tachyon mode with a squared mass in the 
range m\ < m 2 < (and in many cases, it is straightfor- 
ward to prove the existence of a second one in the range 
— oo < m 2 < m 2 ). 

(c) if g + — 0, or if ip' (y + ) = 0, we see that F(Q) = 0: so, 



5 



there is a zero-mode, whose existence is phenomcnolog- 
ically unacceptable, since it mediates (unobserved) long 
range forces, unless the $ wave function is sufficiently 
suppressed on the brane. Precisely, eq. (|3U|) shows that 
^>o {y+) must be bigger than <3>o (0) and than the wave 
function of the massless graviton on the brane. As a con- 
sequence, 3/ would carry long-range forces stronger than 
gravity. 

We insist on the fact that our conclusion doesn't 
rely on any particular assumption concerning the back- 
ground solution. Therefore, any codimension-one com- 
pact model that includes one brane at a fixed point of a 
Z 2 symmetry and a bulk scalar field will have at least a 
tachyon mode, and will be unstable. This agrees with the 
conjecture that was made in Q about the impossibility 
of stabilizing models in which the background behaves at 
some point as in eqs. (|19fl - and proves it in the particular 
case of compact models with only one brane. 



- — c 2 



H(y) 
a (yf 



a (y) 3 dy 



no(») 



k 2 v'o (vY 



lc 2 



a{y) 

f'o (vY 
a(yf 



3 a(y-Y 



H(y- 



a {yf dy 



■ c 2 



(38) 



Now, if g± ^ 0, a zero mode exists if the equations 
IIo (y~) = and Uo(y + ) = can be solved simulta- 
neously with nonvanishing c\ and c 2 . The only case in 
which such condition is verified is if either p' (y-) or 
<Pg (y+) vanish. Therefore, to summarize, zero modes ex- 
ist if and only if at least one of the conditions g + = 0, 
g- =0, p' (?/_) = and ip' (y + ) = holds. In no other 
case massless modes will emerge. 



A condition for stability 



III. TWO BRANES 

In this section we will apply to systems of two branes 
techniques that are similar to the one used in the above 
section. We will deduce some properties of the potentials 
V ((f) and U± (<p) that allow to establish whether any 
tachyon (or zero) modes exist the system. 

To prove the absence of tachyon or zero modes, we 
simply require that the equations 



ff-II m 2 (y_ 



to 2 ty r 



(y- 



o 



(35) 



and 



g + U m 2 (y + ) + to 2 # m 2 (y + ) = (36) 
cannot be solved simultaneously for to 2 < 0. 

A. Zero modes 

The existence of zero modes has been studied in fU . We 
review it here for completeness. First, it is immediate to 
see that if <?+ or g_ vanish, then one of the equations i|35|) 
and (|36[1 with m 2 = is identically verified (whereas the 
other can be used to fix the boundary conditions), indi- 
cating the existence of a zero mode. Suppose now that 
g± 7^ 0. For m 2 = 0, cq. (|llfl can be solved explicitly, 
yielding 



u (y)=c 1 e(y)+c 2 9(y) / d {y)~ 2 dy 



(37) 



Let us now turn to the existence of tachyon modes. 
We will assume that the scale factor a (y) is a mono- 
tonic function of y (therefore the case in which H = 
at some point will not be considered in the present and 
the following subsection), and we choose without loss of 
generality the position of the branes in such a way that 
H (y) < 0. We will also assume that ip' does not van- 
ish at any point (the case in which tp' has zeros will be 
considered in the next subsection). Now suppose that 
g_ > 0. Then cq. I|35|) implies that, for negative m 2 , 
^ m 2 (y_) and H m 2 (y_) must have the same sign. With- 
out loss of generality we can assume it to be positive, 
and suppose {y~) = 1. Then the second of eq. 123|) 
shows that, as long as II m 2 (y) is positive, \t m 2 (jj) > 1. 
But it is easy to see that II m 2 (if) can never get nega- 
tive. Indeed, to get negative, it should cross zero, but for 
n, rl 2 = 0, the first of eqs. IL'-ili shows that (since ip' Q ^ 
by hypothesis) at that point rfll m 2 j dy > 0, so it cannot 
be turning from positive to negative. Thus we conclude 
that for whatever m 2 < 0, ^ m 2 (y + ) and II m 2 (y + ) will 
be positive. As a consequence, if g+ is negative, eq. I|36|) 
cannot be fulfilled for any negative to 2 . Thus, for <?_ > 
and <?+ < 0, no tachyon modes can exist and the system 
turns out to be always stable, as long as the conditions 
H 7^ and ip' ^ hold. From the definition (|17fl we 
see that, if d 2 U±/ dtp 2 — » +oo the system is always stabi- 
lized. In particular, this was the regime analyzed in 0. 
Actually, from the above considerations it follows that 
whatever bulk potential V (p) which does not spoil the 
conditions H ^ and p' ^ can lead to radius stabi- 
lization if d 2 U±/dp 2 ] are positive and large enough. 



that is, 



*o (y) 



Cl 



3 a(yJ] 



H{y-) 



c 2 



H(y) 
a{yf 



C. Instability of other cases 

Keeping the assumption that a is a monotonic function 
of y, and that the position of the branes is chosen in such 



6 



way that H(y) < 0, wc will show here that the three 
criteria formulated in the previous subsection: ip' (y) 7^ 0, 
g- > and < 0, represent a necessary and sufficient 
condition for stability. 

Our strategy will be analogous to the one of the pre- 
vious section, and we will make use of the function 
F (m 2 ) defined in eq. 1|21|) . To begin with, we observe 
that \& m 2 (y_) cannot vanish. Indeed, if this was the 
case the boundary condition (|35[) would impose that also 
^'2 (y~) vanishes (since we know that if g± = a zero 
mode exists, we will assume from now on that g± 7^ 0). 
Then eq. H1()J) would be trivially solved by \E' = 0. We 
will thus be allowed to assume that the function ^ m 2 
appearing in the definition l)21|l is a solution of eq. Hl()(l 
such that \I/ m 2 (y_) = 1, and f ^2 (j/-) is determined by 
eq. G5). 

We can evaluate F (to 2 ) at to 2 = and to 2 — > —00 as 
follows. For m 2 = 0, the solution of eq. (|10|) is given in 
eq. (I38|) . Imposing the boundary condition (|35|) . we get 
exactly the same solution (^o(y),Ilo(y)) as in eas. (l3TJl) . 
and we find that F(Q) is of the sign of g + : 



F(0)=fH|!^ jf- a(yfdy)g + . 



(39) 



Let us then turn to the evaluation of F (m 2 — > —00). In 
the limit of large —to 2 , the equation for v is trivial: 



v" + rn 2 v = . 



(40) 



We can then find the solution for (y) using eq. (|15ll . 
and normalize it to (y_) = 1 and to ea. (|35(l . For 

\m\ (y — y_) 3> 1, the final solution reads 



1 a 3 / 2 (y_) tp' (y) 



g- a 3 / 2 {y) (p' (y-) 

n m 2(y) 



fif- 



Hv-), 

|mj a 3/2 (^_) <^(y) 
3_ a 3 / 2 (?/) ^o(j/_) 



2 / 2 

3# 



'(»-) 



OdTOl" 1 ) 



3 M (y-y-) 
2 ' 

(41) 

e |™| (y-y-) 



Inserting this expression in eq. (|2I|) gives, at leading order 



m to 



lim F (to 



to 



,Jmj a 3/2 (;/-) ^ (y+) el m lfo+-^-> 

<7_ a 3 / 2 (y+) v/ (y_) 2 

(42) 

which is of opposite sign to [g-^g (y+) /Vo (j/-)]- We 
shall now distinguish between various cases: 

(a) if g+ and [g-V'd (y+) /Vo (V-)) nave the same sign, F 
has a zero between m 2 = and to 2 — > —00, so there is 
at least one tachyon mode. 

(b) if g + and [g-(f' (y+) /ip' (y_)] have opposite signs, 



and (/3q vanishes at least in one point j/o m the range 
] U-iV+ii we know that in the limit to 2 — > —00, the so- 
lution for \l/ m 2 (y) must change of sign in the vicinity of 
j/Oj since ip' is in factor of eq. (|41|) . So, we can use the 
same argument as in the single brane case: since (y) 
is strictly positive for m 2 = and goes to negative values 
for to 2 — > —00, there exists a mode to 2 obeying to the 
definition (|32fl . This mode must fulfill condition H34|l . 
and so 



F{m\) = g+n m? (y + ) 



(43) 



is of opposite sign to g+ and to F(0). So, F has a root be- 
tween to 2 = to 2 and to 2 = 0, corresponding to a tachyon 
mode. 

(c) finally, if g + and [g-ip' (y+) /<p' (y_)] have oppo- 
site signs and tp' never vanishes - which implies that 
ip' (y+) /<p' is positive - we can say that: 

• either g + > and g_ < 0, and equation (|41|) shows that 
for any value of y, \I/ m 2 (y) is negative for sufficiently large 
—to 2 . So, we can define the mode to 2 in the same way 
as in case (b), and reach the same conclusion that there 
is at least one tachyon mode. 

• or g + < and g- > 0, and we are exactly in the case 
of the previous subsection, in which we know that there 
is no tachyon mode. 

The cases (a), (b) and (c) include any possible situa- 
tion. So, we conclude that for two-brane systems with a 
monotonic scale factor a(y), the necessary and sufficient 
condition for radion stabilization reads 



9+ < 0, g- > 0, and Vy, <p' (y) / 



(44) 



if the brane positions were chosen in such way that 
H(y) < - the alternative choice H(y) > would lead 
to g+ > 0, g- < 0. 

Criteria that guarantee radius stabilization were stud- 
ied in , using an approach based on a four- 
dimensional effective description of the system. The ex- 
tremization constraint formulated in |ldj | can be rewritten 
in terms of the quantities defined in the present paper as 

(y-f a(y-) v-9- -<p'o (y+f a(y+) y 2 +9+ > o , 

(45) 

where the inequality is saturated when there is a mass- 
less radion. This is indeed a necessary condition, since 
for g_ positive and g + negative the inequality (|45fl al- 
ways holds. On the other hand, the criteria that lead to 
tachyon modes cannot be obtained by the above formula, 
that does not represent, as a consequence, a sufficient 
condition for radius stabilization. Wc attribute this to 
the fact that a four dimensional effective description is 
valid only if the system is actually close to a stable point. 

Before giving an explicit example, we remark that the 
results discussed in the previous section for a single brane 
cannot be obtained as a special case of the present analy- 
sis. It is true that systems with a single brane can be seen 
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as a special case of two brane systems, in two different 
ways. First, one can imagine that one of the two branes 
is fake and has a vanishing tension U- = 0. Second, 
one can impose U- = L/+: then, there is a new symme- 
try with respect to y = (y + + y_)/2, and the system is 
equivalent to a one-brane system plus its mirror image. 
In the first case, however, one would have ip' (y~) = 
and therefore g— would be divergent. In the second case 
one would have that a (y) is not monotonic between the 
two branes. In both cases, the analysis of the present 
section would not hold. This is why we have devoted a 
distinct section to the single brane case. 

D. An example 

As an application of our formalism, let us consider the 
Goldberger-Wise model Q ■ We consider a slight modifi- 
cation of the original model, replacing the quartic terms 
in the brane potentials by quadratic terms, such that the 
bulk and brane potentials read respectively 

U±iv) = + ^ (<P-v ± f . (46) 

In the spirit of the original paper 0] , we have assumed 
for the brane tensions their RS value =f6 / (re 2 t\ . How- 
ever, consistency requires the sum U+ + U— not to be 
positive fLU ]. Therefore, in general either m + or m_ will 
have to be negative. Although negative mass terms are 
present in the five-dimensional lagrangian, this will not 
induce an instability as long as tachyonic modes do not 
appear in the four dimensional effective theory. Alterna- 
tively, we could have chosen to detune the brane tensions, 
thus allowing for a negative value of U+ + U- even for 
positive values of both m+ and m_. 




_1 I 1 1 1 

0.05 0.1 0.15 

M-(y+ - y.) 

FIG. 2: Value of g+ and of the squared mass of the lightest 
Kaluza-Klein mode for the model 1461 . The value of the 
various parameters can be found in the main text. The critical 
value y c below which the system is unstable is such that y c — 
y- ~0.il/x- 1 . 



We now fix £, /j, m_ and u_, solve eqs. (@J in y_ for 
tp and ip', and evolve the bulk equations J5J) until the 
position of the second brane y+ is reached. There we 
read the value of m+ and v+. Therefore, we can deduce 
the behavior of g+ as a function of y+ . In figure [2] such 
behavior is plotted for the following choice of parameters: 
6/(re 2 £ 2 ) = 100, KtL = 1, m_ = 100 //. This implies 
that <7_ ~ 35 fi is positive. The plot shows the presence of 
a critical value y c of y+ such that, for y + < y c , g + has the 
same sign as <?_. Therefore, according the discussion of 
the preceding paragraphs, the system is unstable. On the 
other hand, for y + > y c , the system fulfills the stability 
criteria of the previous section, and no tachyon mode 
should appear. 

As a cross check, in figure |21 we also plot the value 
of the squared mass of the first Kaluza-Klein mode (in 
units of /x 2 ) as a function of y + . Such plot is obtained 
by numerically searching the first zero of the function 
F (m 2 ) defined in JUJ. It is apparent that, when g + 
crosses zero, the squared mass of the first Kaluza-Klein 
mode turns from negative to positive. We thus see that 
the analytical criteria l|44|l arc verified by a full numerical 
computation. Moreover, for the model at hand, we see 
easily that there is a critical distance y c below which the 
model cannot be stable. 



In conclusion, we have seen that the stability and in- 
stability of the system we have considered is essentially 
encoded in the value of g+ and g- (that depend only 
on background quantities). We have shown that mod- 
els including only one brane are unstable, independently 
of the bulk and brane potentials. Zero modes emerge 
if and only if at least one among the quantities </_, <?+, 
tp' (y_) and (p' Q (y+) vanishes. Assuming that H (y) is ev- 
erywhere negative between ?/_ and y+, the condition for 
stability is that Lp' should never vanish in this interval, 
while the quantities g_ and g+ defined in eq. (|17|) should 
be respectively positive and negative. 
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